Partially fluidized shear granular flows: Continuum theory and MD simulations

The continuum theory of partially fluidized shear granular flows is tested and calibrated using two dimensional soft particle molecular dynamics simulations. The theory is based on the relaxational dynamics of the order parameter that describes the transition between static and flowing regimes of granular material. We define the order parameter as a fraction of static contacts among all contacts between particles. We also propose and verify by direct simulations the constitutive relation based on the splitting of the shear stress tensor into a``fluid part'' proportional to the strain rate tensor, and a remaining ``solid part''. The ratio of these two parts is a function of the order parameter. The rheology of the fluid component agrees well with the kinetic theory of granular fluids even in the dense regime. Based on the hysteretic bifurcation diagram for a thin shear granular layer obtained in simulations, we construct the ``free energy'' for the order parameter. The theory calibrated using numerical experiments with the thin granular layer is applied to the surface-driven stationary two dimensional granular flows in a thick granular layer under gravity.Comment: 20 pages, 19 figures, submitted to Phys. Rev.

Vector lattice model for stresses in granular materials

A vector lattice model for stresses in granular materials is proposed. A two dimensional pile built by pouring from a point is constructed numerically according to this model. Remarkably, the pile violates the Mohr Coulomb stability criterion for granular matter, probably because of the inherent anisotropy of such poured piles. The numerical results are also compared to the earlier continuum FPA model and the (scalar) lattice $q$-model

Granular packings with moving side walls

The effects of movement of the side walls of a confined granular packing are studied by discrete element, molecular dynamics simulations. The dynamical evolution of the stress is studied as a function of wall movement both in the direction of gravity as well as opposite to it. For all wall velocities explored, the stress in the final state of the system after wall movement is fundamentally different from the original state obtained by pouring particles into the container and letting them settle under the influence of gravity. The original packing possesses a hydrostatic-like region at the top of the container which crosses over to a depth-independent stress. As the walls are moved in the direction opposite to gravity, the saturation stress first reaches a minimum value independent of the wall velocity, then increases to a steady-state value dependent on the wall-velocity. After wall movement ceases and the packing reaches equilibrium, the stress profile fits the classic Janssen form for high wall velocities, while it has some deviations for low wall velocities. The wall movement greatly increases the number of particle-wall and particle-particle forces at the Coulomb criterion. Varying the wall velocity has only small effects on the particle structure of the final packing so long as the walls travel a similar distance.Comment: 11 pages, 10 figures, some figures in colo

Confined granular packings: structure, stress, and forces

The structure and stresses of static granular packs in cylindrical containers are studied using large-scale discrete element molecular dynamics simulations in three dimensions. We generate packings by both pouring and sedimentation and examine how the final state depends on the method of construction. The vertical stress becomes depth-independent for deep piles and we compare these stress depth-profiles to the classical Janssen theory. The majority of the tangential forces for particle-wall contacts are found to be close to the Coulomb failure criterion, in agreement with the theory of Janssen, while particle-particle contacts in the bulk are far from the Coulomb criterion. In addition, we show that a linear hydrostatic-like region at the top of the packings unexplained by the Janssen theory arises because most of the particle-wall tangential forces in this region are far from the Coulomb yield criterion. The distributions of particle-particle and particle-wall contact forces $P(f)$ exhibit exponential-like decay at large forces in agreement with previous studies.Comment: 11 pages, 11 figures, submitted to PRE (v2) added new references, fixed typo

Stress Transmission through Three-Dimensional Ordered Granular Arrays

We measure the local contact forces at both the top and bottom boundaries of three-dimensional face-centered-cubic and hexagonal-close-packed granular crystals in response to an external force applied to a small area at the top surface. Depending on the crystal structure, we find markedly different results which can be understood in terms of force balance considerations in the specific geometry of the crystal. Small amounts of disorder are found to create additional structure at both the top and bottom surfaces.Comment: 9 pages including 9 figures (many in color) submitted to PR

Effect of boundaries on the force distributions in granular media

The effect of boundaries on the force distributions in granular media is illustrated by simulations of 2D packings of frictionless, Hertzian spheres. To elucidate discrepancies between experimental observations and theoretical predictions, we distinguish between the weight distribution {\cal P} (w) measured in experiments and analyzed in the q-model, and the distribution of interparticle forces P(f). The latter one is robust, while {\cal P}(w) can be obtained once the local packing geometry and P(f) are known. By manipulating the (boundary) geometry, we show that {\cal P}(w) can be varied drastically.Comment: 4 pages, 4 figure

Modeling relaxation and jamming in granular media

We introduce a stochastic microscopic model to investigate the jamming and reorganization of grains induced by an object moving through a granular medium. The model reproduces the experimentally observed periodic sawtooth fluctuations in the jamming force and predicts the period and the power spectrum in terms of the controllable physical parameters. It also predicts that the avalanche sizes, defined as the number of displaced grains during a single advance of the object, follow a power-law, $P(s)\sim s^{-\tau}$, where the exponent is independent of the physical parameters

The effects of grain shape and frustration in a granular column near jamming

We investigate the full phase diagram of a column of grains near jamming, as a function of varying levels of frustration. Frustration is modelled by the effect of two opposing fields on a grain, due respectively to grains above and below it. The resulting four dynamical regimes (ballistic, logarithmic, activated and glassy) are characterised by means of the jamming time of zero-temperature dynamics, and of the statistics of attractors reached by the latter. Shape effects are most pronounced in the cases of strong and weak frustration, and essentially disappear around a mean-field point.Comment: 17 pages, 19 figure

Statistics of the contact network in frictional and frictionless granular packings

Simulated granular packings with different particle friction coefficient mu are examined. The distribution of the particle-particle and particle-wall normal and tangential contact forces P(f) are computed and compared with existing experimental data. Here f equivalent to F/F-bar is the contact force F normalized by the average value F-bar. P(f) exhibits exponential-like decay at large forces, a plateau/peak near f = 1, with additional features at forces smaller than the average that depend on mu. Computations of the force-force spatial distribution function and the contact point radial distribution function indicate that correlations between forces are only weakly dependent on friction and decay rapidly beyond approximately three particle diameters. Distributions of the particle-particle contact angles show that the contact network is not isotropic and only weakly dependent on friction. High force-bearing structures, or force chains, do not play a dominant role in these three dimensional, unloaded packings.Comment: 11 pages, 13 figures, submitted to PR

Properties of layer-by-layer vector stochastic models of force fluctuations in granular materials

We attempt to describe the stress distributions of granular packings using lattice-based layer-by-layer stochastic models that satisfy the constraints of force and torque balance and non-tensile forces at each site. The inherent asymmetry in the layer-by-layer approach appears to lead to an asymmetric force distribution, in disagreement with both experiments and general symmetry considerations. The vertical force component probability distribution is robust and in agreement with predictions of the scalar q model while the distribution of horizontal force components is qualitatively different and depends on the details of implementation.Comment: 18 pages, 12 figures (with subfigures), 1 table. Uses revtex, epsfig,subfigure, and cite. Submitted to PRE. Plots have been bitmapped. High-resolution version is available. Email [email protected] or download from http://rainbow.uchicago.edu/~mbnguyen/research/vm.htm